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Abstract 

A critical assessment is given of the theoretical uncertainty in the predicted cross-sections 
for large-angle Bhabha scattering at LEP 1, with or without t-channel subtraction. To 
this end a detailed comparison is presented of the results obtained with the programs 
ALIBABA and TOPAZO. Differences in the implementation of the radiative corrections 
and the effect of missing higher-order terms are critically discussed. 
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1 Introduction 



In 1995 the Z phase of LEP has come to an end and at present the ultimate analysis of the data is 
imminent. This involves in particular the completion of the line-shape analysis, including the final 
LEP energy calibration. Consequently, the safest possible estimate of the theoretical accuracy is of 
the upmost importance. It should be noted that the LEP 1 data (1990-1995) have been taken in 
the energy ) range \^/s — M^\ < 3 GeV and consist of the hadronic and leptonic cross-sections, 
the leptonic forward-backward asymmetries, the various polarization asymmetries, the partial widths, 
and the quark forward-backward asymmetries. All this makes it mandatory to assess the theoretical 
precision of the available programs for different channels and for energies up to ~ 3 GeV away from 
the resonance. 

In this note we focus on the electron (Bhabha) channel. Bhabha scattering is measured with 
remarkable precision at LEP 1/SLC in two complementary kinematical regions: small and large scat- 
tering angles. The former plays a crucial role in determining the luminosity. The latter is essential for 
extracting the Z-boson properties. The main message of this note concerns an update of the theoret- 
ical precision in the large-angle regime. This precision depends on the beam energy and on the event 
selection used, the worst case being at a few GeV above the Z resonance. There the different programs 
are found to deviate by as much as 1%. This should be contrasted with the expected experimental 
systematic errors |j|] , displayed in Table || (an example of the statistical errors is given in Table . 
In this note we present a detailed analysis of the observed deviations and translate this into estimates 
for the theoretical precision. 





ALEPH 


DELPHI 


L3 


OPAL 




0.30/0.30/0.31% 
0.0018/0.0019/0.0020 


0.59/0.54/0.75% 
0.0025/0.0022/0.0021 


0.30/0.23/1.0% 
0.003/0.003/0.01 


0.23/0.24/ - % 
0.0016/0.0016/0.002 



Table 1: The (preliminary) experimental systematic errors for the years 1993/1994/1995 at the Z 
peak, not including the common uncertainty due to the LEP energy calibration. 



2 Comparison of ALIBABA and TOPAZO 

In order to assess the theoretical uncertainties in the predictions for large-angle Bhabha scattering 
(LABS), we first perform a detailed numerical comparison of the programs ALIBABA |^] and TOPAZO 
[^. Previous comparisons can be found in |5| and also in [^]. The input parameters chosen for the 
comparison are: = 91.1863 GeV, rrit = 175.6 GeV, = 300 GeV, and as(M|) = 0.118. 

The other Standard Model parameters we take from Ref. 0]. The acceptance cuts for the final-state 
particles consist of a minimum energy for both particles (1 GeV), an angular acceptance for the 
electron (40° < < 140°), and a maximum acollinearity angle (10° or 25°). For the energy we take 
the characteristic LEP 1 energies: ^/s = 88.45, 89.45, 90.20, 91.19, 91.30, 91.95, 93.00, and 93.70 GeV. 

In Table ^ we present the comparison for the cross-sections (a) and forward-backward asymmetries 
(Aps)- Besides the results for the full Bhabha process, we also show the pure s-channel contributions. 
The reason for that is twofold. First, the s-channel branch of TOPAZO allows a systematic inclusion of 
the full 0(0^) QED corrections [^, which are only present in leading- log approximation in ALIBABA. 
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1990 




[GeV] 


88.223 


91.215 


94.202 


results 


stat 


error 


3.8% 


2.5% 


3.1% 


1991 




[GeV] 


88.464 


91.207 


93.701 


results 


stat 


error 


3.3% 


2.2% 


3.1% 


1992 




[GeV] 




91.270 




results (SICAL) 


stat 


error 




1.6% 




1993 


V~s 


[GeV] 


89.43204 


91.18718 


93.01527 


results 


stat 


error 


1.2% 


1.6% 


1.1% 


1994 




[GeV] 




91.19677 




results 


stat 


error 




0.7% 




1995 




[GeV] 


89.43952 


91.28190 


92.96812 


results 


stat 


error 


1.2% 


2.2% 


1.0% 



Table 2: Sample of statistical errors for the electron cross-section from ALEPH. For the various runs 
the highest and lowest energies are shown, as well as the energy closest to the Z peak. 

In this way some of the statements that will be made for the full Bhabha process can be checked in 
the s-channel mode. Secondly, the so-called i-channel subtraction procedure for extracting the pure 
s-channel part is still very popular among the experimental Collaborations. As such comparisons of 
all components of LABS are well deserving the effort. 

The numerical precision of the TOPAZO numbers is 0.00001 in ApB and 0.001% (relative precision) 
in a. For ALIBABA the numerical precision is much less in view of the five-dimensional VEGAS 
integration, i.e. 0.0006 in Aps and 0.05% in a. The technical precision of the results, however, will 
not affect our conclusions. Whereas the observed LABS differences between TOPAZO and ALIBABA 
are largely acceptable for a maximum acollinearity angle of 25°, this is certainly not true for 10°. In 
that case large deviations up to 1% show up for energies above the Z resonance, as already discussed 
ini. 

In Table § we try to quantify where this difference is coming from. In this context it is important 
to note that the main conceptual difference between TOPAZO and ALIBABA in the Bhabha channel 
lies in the implementation of the non-leading-log QED corrections. Both programs are based on the 
structure- function method for calculating the (dominant) leading- log corrections [oc (oL/tt)", with 
n = 1,2 and L = log(,s/r?T,g)]. When it comes to the non-leading-log corrections both programs use 
quite different approaches. TOPAZO was designed to be an efficient fitter of realistic observables 
around the Z resonance, where LABS is effectively dominated by the annihilation channel. The 
structure functions in TOPAZO are based on the iterative solution of the master evolution equation, 
which accounts for the well-known second-order electron form factor, and therefore they reproduce 
the sub-leading terms [oc (a/7r)"L'", with n = 1,2 and n > m] for the s channel.^] This procedure, 
however, does not reproduce the (unknown) correct answer for the sub-leading terms in the t channel. 

In TOPAZO a part of the sub-leading terms in the t channel enters through a final-state-radiation 
factor, which is subsequently convoluted with the structure functions. In this way part of the 0{a) 

^In a realistic (i.e. not fully-extrapolated) set-up some 0{a) hard-photonic non-log corrections are missing in the s 
channel. However, their effect is tiny around the Z resonance 
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LEP 1 energy in GeV 


88.45 


89.45 


90.20 


91.19 


91.30 


91.95 


93.00 


93.70 


maximum acollinearity angle: 10° 


0-''+* (T) 
(A) 


457.08 
457.71 
-0.14% 


644.86 
644.78 
+0.01 % 


912.06 
911.43 
+0.07% 


1185.70 
1184.59 
+0.09% 


1164.82 
1163.71 
+0.10% 


873.50 
876.40 
-0.33% 


476.64 
480.23 
-0.75 % 


351.80 
355.31 
-1.00% 


(T) 

(A) 
(T-A) 


+0.4448 
+0.4454 
-0.0006 


+0.3411 
+0.3409 
+0.0002 


+0.2492 
+0.2489 
+0.0003 


+0.1386 
+0.1389 
-0.0003 


+0.1298 
+0.1301 
-0.0003 


+0.1008 
+0.1020 
-0.0012 


+0.1298 
+0.1315 
-0.0017 


+0.1788 
+0.1818 
-0.0030 


(T) 
(A) 
(S) 


172.94 
173.60 
-0.38% 


331.55 
332.09 
-0.16% 


590.93 
590.72 
+0.04% 


994.27 
991.93 
+0.24% 


998.32 
996.78 
+0.15% 


820.80 
821.13 
-0.04% 


461.49 
463.35 
-0.40 % 


329.61 
331.49 
-0.57% 


A'fb (T) 
(A) 
(T-A) 


-0.2202 
-0.2209 
+0.0007 


-0.1380 
-0.1386 
+0.0006 


-0.0761 
-0.0774 
+0.0013 


+0.0004 
-0.0008 
+0.0012 


+0.0081 
+0.0072 
+0.0009 


+0.0487 
+0.0485 
+0.0002 


+0.0980 
+0.0977 
+0.0003 


+0.1225 
+0.1225 
+0.0000 


maximum acollinearity angle: 25° 


a'+' (T) 
(A) 


485.17 
484.05 
+0.23% 


674.89 
673.91 
+0.15% 


945.00 
944.73 
+0.03 % 


1221.13 
1220.49 
+0.05 % 


1200.16 
1199.51 
+0.05 % 


905.25 
907.15 
-0.21% 


503.79 
504.73 
-0.19% 


377.59 
378.40 
-0.21% 


Ap-^ (T) 

(A) 
(T-A) 


+0.4605 
+0.4576 
+0.0029 


+0.3554 
+0.3521 
+0.0033 


+0.2613 
+0.2596 
+0.0017 


+0.1501 
+0.1484 
+0.0017 


+0.1417 
+0.1402 
+0.0015 


+0.1175 
+0.1173 
+0.0002 


+0.1584 
+0.1580 
+0.0004 


+0.2136 
+0.2130 
+0.0006 


(T) 
(A) 
(S) 


176.31 
177.43 
-0.64% 


336.84 
338.44 
-0.48 % 


599.25 
601.29 
-0.34% 


1007.03 
1008.13 
-0.11% 


1011.10 
1011.67 
-0.06 % 


831.43 
833.29 
-0.22% 


468.16 
469.57 
-0.30 % 


334.94 
335.92 
-0.29% 


A'fb (T) 
(A) 
(T-A) 


-0.2235 
-0.2227 
-0.0008 


-0.1404 
-0.1406 
+0.0002 


-0.0777 
-0.0777 
+0.0000 


-0.0004 
-0.0007 
+0.0003 


+0.0073 
+0.0069 
+0.0004 


+0.0480 
+0.0481 
-0.0001 


+0.0967 
+0.0976 
-0.0009 


+0.1203 
+0.1210 
-0.0007 



Table 3: Comparison of TOPAZO (T) and ALIBABA (A) for the cross-section (in pb) and the forward-^ 
backward asymmetry. The TOPAZO results do not include initial-state pair production. The full 
Bhabha results are indicated by the superscript "s + t", the s-channel contributions by "s". The 
quantity 5 stands for the relative deviation 100% . (T — A)/T. The input parameters can be found in 
the text. 

non-log terms are implemented, as well as a subset of sub-leading second-order terms oc a L/tt^. These 
sub-leading higher-order terms are not present in ALIBABA. However, in the ALIBABA approach 
the full set of non-log 0{a) corrections are determined by means of the five-dimensional VEGAS 
integration, needed for handling the radiative process. These non-log corrections are listed in the 
NLa entries in Table ^. The entries A^l in Table ^ quantify the deviations caused by considering 
non-leading terms only from final-state radiation. In conclusion, the main difference in the treatment 
of the QED corrections in both programs can be estimated by Cp {a^~^^ — NL/^) — NL^. Here the 
0{a) final-state-radiation factor Cp is given by —0.0165 (—0.0052) for a maximum acollinearity angle 
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LEP 1 energy in GeV 


88.45 


89.45 


90.20 


91.19 


91.30 


91.95 


93.00 


93.70 


maximum acollinearity angle: 10° 


NLa [pb] 


-8.05 


-11.54 


-15.45 


-17.80 


-17.52 


-10.48 


-4.71 


-2.88 


Afsr 
S - ApsR 


-0.33 % 
+0.08 % 
-0.22% 


-0.39% 
+0.11% 
-0.10% 


-0.54% 
+0.02 % 
+0.05 % 


-0.64% 
-0.17% 
+0.27% 


-0.60% 
-0.17% 
+0.26% 


-0.62 % 
-0.48 % 
+0.14% 


-0.42% 
-0.69% 
-0.06% 


-0.44% 
-0.86 % 
-0.14% 


maximum acollinearity angle: 25° 


NLa [pb] 


-4.14 


-4.59 


-4.17 


-3.74 


-3.54 


-1.48 


-1.89 


-1.49 


Anl 

ApSR 

6 - ApsR 


+0.23 % 
+0.34% 
-0.11% 


+0.02 % 
+0.16% 
-0.01% 


-0.24% 
-0.08 % 
+0.11% 


-0.35% 
-0.22% 
+0.27% 


-0.35% 
-0.23% 
+0.28% 


-0.39% 
-0.36 % 
+0.15% 


-0.05% 
-0.15% 
-0.04% 


+0.02 % 
-0.13% 
-0.09 % 



Table 4: An analysis of the differences between TOPAZO and ALIBABA for large-angle Bhabha 
scattering. The various entries are explained in the text. 

of 10° (25°). In Table § we present this estimate in the form of a relative correction factor (Apsr) 
with respect to the full Bhabha cross-section. Comparing the estimate with the actual deviation 6 
in Table ^, a systematic shift {S — Apsr) is observed that is roughly independent of the maximum 
acollinearity angle. The resulting shift can be attributed to the fact that the ALIBABA weak library 
is not up to date. This is supported by a direct comparison with the Bhabha branch of the program 
ZFITTER Q . Another thing to note is that for a maximum acollinearity angle of 10° both TOPAZO 
and ALIBABA seem to miss terms of the order of 0.3-0.6% (cf. A^l and A^l — ApsR in Table §). 
Below and near the Z resonance both effects go in the same direction, compensating each other in 
the difference. This leads to a somewhat misleading agreement between the programs. Above the 
resonance the effects have different signs, leading to an enhancement of the deviations. The relative 
deviation in the non s-channel components {a^~^* — a'^), which is practically negligible below the 
resonance, grows up to ~ 10% on the high-energy side. There, however, the non s-channel cross- 
sections are rapidly decreasing (e.g. roughly 23 pb at -y/i = 93.70 GeV, i.e. 6.5% of full LABS, versus 
284 pb at ^ = 88.45 GeV, i.e. 62% of full LABS). 

3 Theoretical error estimates for large-angle Bhabha scattering 

Having established the main sources of the differences between TOPAZO and ALIBABA, we can now 
address the question of the error in the theoretical predictions for large-angle Bhabha scattering. 
In the experimental analyses either the full Bhabha cross-section is used or merely the s-channel 
contributions. The latter are obtained through t-channel subtraction, i.e. by subtracting the non 
s-channel contributions that involve i-channel gauge-boson exchange. The subtraction procedure is 
aimed at reducing the full LABS to a simple annihilation process to be subsequently analyzed by some 
up-to-date (from the point of view of the WEAK/QCD library) program. 

Correspondingly we will present the theoretical errors for both procedures. The bulk of the the- 
oretical errors are due to missing QED corrections: missing non- log 0{a) corrections in TOPAZO, 
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missing higher-order sub- leading effects from the structure- function convolution of the non-log 0{a) 
corrections in ALIBABA, missing effects in both programs from the (unknown) sub-leading terms 
in the structure functions, and missing initial-state pair-production effects in ALIBABA. In addition 
there is the usual calorimetric measurement problem, which arises if the cross-section is not inclusive 
in the energy of the outgoing fermions. While with the present cuts the effect is estimated to be 
negligible by TOPAZO 0.1%), for high-energy thresholds the contribution will grow considerably. 

In general one should remember that an absolute test of the QED corrections does not exist. QED 
corrections are convoluted with - in principle - different kernel cross-sections, so that the shift in the 
absolute QED corrections also depends on the differences in the non-QED parts. 

Before coming to all these QED effects, we first estimate the uncertainty in the weak sector. From 
a detailed study with TOPAZO it follows that different weak options influence the cross-sections by at 
most 0.06% and the forward-backward asymmetry by at most 0.0002. In this respect the quality of the 
estimate has not changed since the work done in ||5|, since the upgrading in TOPAZO and ZFITTER 
have been constantly cross-checked. The uncertainty in ALIBABA roughly equals the ApsR ~ ^ 
shifts displayed in Table ^ varying between 60 and 1000 GeV, and M^, nit, and asiM"^) within 
one experimental standard deviation, leads to variations of the cross-sections and asymmetries that 
closely resemble the ones observed with TOPAZO. The ALIBABA variations have the same sign as 
the TOPAZO ones and do not differ by more than 0.06% for the cross-section and 0.0006 for the 
forward-backward asymmetry, which is of the same size as the numerical error in the ALIBABA 
results. 

In the context of the QED errors, the ALIBABA error owing to the absence of initial-state pair- 
production (ISPP) effects can be estimated in a straightforward way with the help of TOPAZO. For 
the eight LEP 1 energy points we flnd for the full LABS cross-section: —0.23%, —0.25%, —0.26%, 
-0.25%, -0.24%, -0.19%, -0.04%, +0.07%. Here the final-state e+e" pair is required to have an 
invariant mass that exceeds 50% of the total energy ^/s (i.e. Zmin = 0.25). It should be noted that 
among all radiative corrections the ISPP appears as the most questionable. While relatively safe 
around the resonance, different implementations of ISPP start to register some disagreement at higher 
energies |10| and strongly depend on kinematical cuts, in other words on the precise separation between 
2-fermion and 4-fermion physics. 

The missing sub-leading higher-order terms from the convolution of the non-log 0{a.) corrections 
(MSLl) in ALIBABA is harder to estimate. Here we deploy two methods. As can be read off from the 
NLyi entries in Table ^, the non-log 0{a) corrections in ALIBABA appear in two different shapes. For 
a maximum acollinearity angle of 10° the NL^ contribution exhibits a resonance-like behaviour, with 
a pronounced peak around ^/s = M^. In that case the leading-log corrections to this contribution 
should closely resemble the leading-log corrections to the lowest-order cross-section. A natural way 
of estimating these higher-order effects is by dividing the NL^i contributions by the lowest-order 
cross-sections (in 'dressed-Born' form), taking the maximum value of this ratio, and subsequently 
multiplying the leading-log corrections in ALIBABA by this number. Note that in principle the sign 
of the missing terms will be fixed. For a maximum acollinearity angle of 25° the NL^ contribution is 
rather flat, with a sudden jump just above the Z resonance. In that case the best way to estimate the 
higher-order effects is by simply multiplying the maximum value for |NLyi| by the factor (4aL/7r). Here 
an additional factor of two is added for safety. The sign of the missing terms is only determined for 
the three highest energy points, i.e. the ones after the jump. Since the structure-function convolution 
probes the variation in the NL^^ contributions at lower energies, a negative sign is expected for the 
higher-order corrections in these three energy points. Note that this second method also applies to 
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the non s-channel cross-sections needed for the t-channel subtraction, irrespective of the maximum 
acollinearity angle. Finally, the effect from the (unknown) sub-leading terms in the structure functions 
(MSL2) can be estimated by multiplying the leading-log corrections by the factor (4a/7r). Again a 
factor of two is added for safety. Note that the sign of the MSL2 terms is not fixed. 

In Table ^ we summarize the estimates for the theoretical errors of the ALIBABA program.^ In 
the total error the various estimates are added in quadrature. For the unsubtracted cross-section we 





LEP 1 energy in GeV 




88.45 


89.45 


90.20 


91.19 


91.30 


91.95 


93.00 


93.70 






maximum acollinearity angle: 


10° 






weak 


0.17% 


0.06 % 


0.08 % 


0.27% 


0.27% 


0.15% 


0.05 % 


0.12% 




0.57 pb 


0.63 pb 


0.64 pb 


0.39 pb 


0.33 pb 


0.33 pb 


0.33 pb 


0.33 pb 


pairs 


0.23 % 


0.25 % 


0.26 % 


0.25 % 


0.24% 


0.19% 


0.04% 


0.07% 




0.65 pb 


0.78 pb 


0.83 pb 


0.48 pb 


0.40 pb 


0.32 pb 


0.07pb 


0.12 pb 


MSLl 


0.34 % 


0.42 % 


0.47% 


0.39 % 


0.36 % 


0.13% 


0.22 % 


0.35 % 




1.36 pb 


1.36 pb 


1.36 pb 


1.36 pb 


1.36 pb 


1.36 pb 


1.36 pb 


1.36 pb 


MSL2 


0.23 % 


0.28 % 


0.32 % 


0.26 % 


0.24% 


0.24% 


0.24% 


0.24% 




0.46 pb 


0.54 pb 


0.50 pb 


0.31 pb 


0.42 pb 


0.76 pb 


0.52 pb 


0.35 pb 


total 


0.5% 


0.6% 


0.6% 


0.6% 


0.6% 


0.4% 


0.3% 


0.4% 




1.7pb 


1.8 pb 


1.8 pb 


1.5 pb 


1.5 pb 


1.6 pb 


1.5 pb 


1.4 pb 






maximum acollinearity angle: 


25° 






weak 


0.17% 


0.06 % 


0.08 % 


0.27% 


0.27% 


0.15% 


0.05 % 


0.12% 




0.61 pb 


0.67 pb 


0.69 pb 


0.42 pb 


0.38 pb 


0.38 pb 


0.38 pb 


0.38 pb 


pairs 


0.23 % 


0.25 % 


0.26 % 


0.25 % 


0.24% 


0.19% 


0.04 % 


0.07% 




0.71 pb 


0.84 pb 


0.89 pb 


0.53 pb 


0.45 pb 


0.36 pb 


0.08 pb 


0.13 pb 


MSLl 


0.21 % 


0.15% 


0.11% 


0.08 % 


0.09 % 


0.11% 


0.20 % 


0.27% 




0.97 pb 


0.97 pb 


0.97pb 


0.97 pb 


0.97 pb 


0.97pb 


0.97pb 


0.97 pb 


MSL2 


0.17% 


0.24% 


0.28 % 


0.24 % 


0.22% 


0.22 % 


0.22 % 


0.27% 




0.26 pb 


0.35 pb 


0.30 pb 


0.50 pb 


0.62 pb 


0.95 pb 


0.71 pb 


0.53 pb 


total 


0.4% 


0.4% 


0.4% 


0.4% 


0.4% 


0.3% 


0.3% 


0.4% 




1.4 pb 


1.5 pb 


1.5 pb 


1.3 pb 


1.3 pb 


1.5 pb 


1.3 pb 


1.2 pb 



Table 5: Error estimates for the large-angle Bhabha cross-section as predicted with ALIBABA. The 
numbers in the first row of every entry correspond to the unsubtracted cross-sections and are given 
relative to o"*"*"* (see Table ^). The numbers in the second row (given in pb) correspond to the non 
s-channel components needed for the i-channel subtraction. 

take the average of the two 5 — ApsR shifts, displayed in Table ^, as a measure of the error in the 
weak corrections. For the non s-channel cross-section a constant uncertainty of 0.2% is assumed. At 
this point one should bear in mind that the non s-channel cross-sections exhibit a strong cancellation 

^We have also derived the non s-channel errors for the ALEPH angular acceptance, involving a 20° maximum 
acollinearity angle and 45° < ^(e~,ebcam) < 155°. The errors read: 1.1 pb, 1.2 pb, 1.2 pb, 1.1 pb, 1.1 pb, 1.2 pb, 1.1 pb, 
and 1.0 pb. These errors are smaller than the ones in Table ^ since also the non s-channel cross-sections themselves are 
smaller. The errors on the full LABS cross-section are given by; 0.5%, 0.4%, 0.4%, 0.5%, 0.5%, 0.4%, 0.4%, and 0.5%. 
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between the t-channel and s-t interference contributions at the three highest energy points (see Ta- 
ble ^). In order to avoid underestimating the weak and ISPP errors at these energies, we therefore use 
the cross-section at 91.30 GeV in the error estimates. The error estimates for the non s-channel cross- 
sections are given in pb, rather than relative to the full fj*"*"* — results. In this way the errors are 
independent of the amount of cancellation between the t-channel and s-t interference contributions. 
Note that the non s-channel errors are not derived from the errors for the full Bhabha cross-section. 
They are derived directly from a^~^^ — u*, eliminating in this way the correlation between the errors 
in (7*+* and a'^. 

As a cross-check of the above error-estimate procedure, we have also derived the pure s-channel 
errors. For a maximum acollinearity angle of 10° we find 0.4%, 0.4%, 0.4%, 0.5%, 0.5%, 0.3%, 0.4%, 
and 0.5%; for 25° the errors are 0.7%, 0.6%, 0.5%, 0.3%, 0.3%, 0.4%, 0.4%, and 0.4%. This is in good 
agreement with the observed s-channel deviations in Table ^. 

The main approximation for LABS in TOPAZO is the leading-logarithmic one in the non s-channel 
components. The QED theoretical error in TOPAZO is therefore growing when enlarging the angular 
acceptance to smaller angles, or in general in all situations where the non s-channel terms are increased. 

In order to give an estimate for the theoretical error in TOPAZO we proceed in the following - 
conservative - way. First of all, we assign an overall 0.06% error to the weak corrections. This is 
not the consequence of a mistreatment of the corrections, but a safe upper bound on the uncertainty 
arising from different implementation schemes. 

Next, the ISPP as implemented in TOPAZO is strictly speaking only valid for s-channel processes. 
If Ap is the relative ISPP effect, then Apfi* is the correct s-channel contribution. The missing terms 
are proportional to the pure t-channel and s~t interference components. However, we refrain from 
assigning to LABS an overall uncertainty of =bAp (u^^* — o"*), since the difference 0"*+* — becomes 
very small above the peak, while the individual t-channel and s-t interference terms can be as large 
as below the peak. Clearly this cancellation cannot be transferred to the error estimates. We use 
two different procedures. First we proceed by assigning to the cross-section an error zbApUmax, where 
<7max = maxjcr*, cj'*"^* — a^} . The obtained error is roughly half of Apcr'*"'"* below the resonance and 
almost coincides with ApO"*"^* above it. Alternatively we compute dabs) the sum of the t-channel term 
and the absolute value of the s-t interference, and assign an error ztApCabs- In the two procedures 
the errors read: 0.14%, 0.13%, 0.17%, 0.21%, 0.21%, 0.18%, 0.04%, 0.07% and 0.17%, 0.14%, 0.11%, 
0.05%, 0.05%, 0.08%, 0.03%, 0.07%. There are appreciable differences only around the peak. In the 
summary we will report the average of the two procedures. 

Similarly the MSL2 effect in TOPAZO is only active for the t-channel and s-t interference con- 
tributions. For the associated uncertainty we take again the average of MSL2(Tmax and MSL2(Tabs) 
with MSL2 as estimated by ALIBABA. We have also performed a consistency check by computing 
the TOPAZO non s-channel cross-sections with the change L — > (lib4^)L in the definition of /3, 
appearing in the structure functions. The resulting effect is roughly equal to the one obtained in the 
(Tabs method. Nevertheless we prefer to be more conservative and to average with the error from the 
(Tmax method. 

Another source of theoretical error remains in the missing non- log 0{a) corrections (MNL) from 
QED radiation. TOPAZO has some of these non-log terms, which have been included into the structure 
functions through the exact second-order s-channel vertex. We have computed this contribution, Aj^j^ , 
and found that it is largely due to the 7r^/3 — 2 term in the K-factor. Next we assume that — AtvLj as 
taken from ALIBABA, is the exact contribution and we estimate the uncertainty to be — Ajvl — A]^^. 
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In Table | we summarize the estimates for the theoretical errors of the TOPAZO program. In the 
total error the various estimates are again added in quadrature. 





LEP 1 energy in GeV 


88.45 


89.45 


90.20 


91.19 


91.30 


91.95 


93.00 


93.70 


maximum acollinearity angle: 10° 


weak 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


pairs 


0.16% 


0.14% 


0.14% 


0.13% 


0.13% 


0.13% 


0.04% 


0.07% 


MSL2 


0.16% 


0.15% 


0.17% 


0.14% 


0.13% 


0.17% 


0.21 % 


0.23 % 


MNL 


0.00 % 


0.07% 


0.23 % 


0.33 % 


0.29 % 


0.28 % 


0.03 % 


0.01 % 


total 


0.2% 


0.2% 


0.3% 


0.4% 


0.4% 


0.4% 


0.2% 


0.3% 


maximum acollinearity angle: 25° 


weak 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


0.06 % 


pairs 


0.17% 


0.14% 


0.14% 


0.14% 


0.13% 


0.13% 


0.04% 


0.07% 


MSL2 


0.12% 


0.13% 


0.15% 


0.13% 


0.12% 


0.15% 


0.19% 


0.25 % 


MNL 


0.65 % 


0.41 % 


0.12% 


0.00 % 


0.00 % 


0.01 % 


0.42 % 


0.56 % 


total 


0.7% 


0.5% 


0.3% 


0.2% 


0.2% 


0.2% 


0.5% 


0.6% 



Table 6: Error estimates for the large-angle Bhabha cross-section as predicted with TOPAZO. 

Finally we come to the theoretical errors on the forward-backward asymmetry. As can be seen 
from Table ^ the forward-backward asymmetry shows the typical t-channel effect, which makes 
it unique among leptonic asymmetries. At the peak we observe a T— A absolute difference of —0.0003 
(-F0.0017) for an average asymmetry of -h0.1388 (-hO.1493) at ^acoU = 10° (25°). This deteriorates up 
to a -0.0030 (-F0.0033) difference for an average asymmetry of -^0.1803 (-hO.3538) at ^ = 93.70 GeV 
(89.45 GeV) and 

^acoii — 10° (25°). The s-channel asymmetry, -A^p^, is instead very small in the Z 
peak region. For this quantity the disagreement is globally contained within an absolute difference 
of 0.0013. This is no surprise, since it has been shown in Q that differences among programs are 
weak-dominated around the peak and QED-dominated far from the peak only. For Ap'^'^ we assign 
to both codes an error ± A^^ , with A^^ given by the half-difference of the predictions (i.e. half of the 
T-A row of table |). 

4 Conclusions 

The main emphasis of this note was on presenting the safest possible estimate for the theoretical 
accuracy of the large-angle Bhabha scattering calculations. As usual, in estimating theoretical uncer- 
tainties one proceeds by comparing the results of different programs starting from a common set of 
input parameters and kinematical cuts. In our case we have used the predictions of ALIBABA and 
TOPAZO. The registered differences give a rough idea of the uncertainty associated with different 
implementations of radiative corrections. Successively one tries to estimate the internal accuracy of 
the programs by deriving bounds on the effects of missing parts of the calculations. 

Our analysis shows a substantial agreement between the two programs for larger acollinearity 
cuts, with a deterioration towards small acollinearity cuts and energies on the high-energy side of the 
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resonance. For a maximum acollinearity angle of 10° we find deviations of —0.14%, +0.09%, —1.00% 
at - 2.5 GeV, M^, +2.5 GeV for the full Bhabha cross-section (-0.38%, +0.24%, -0.57% for 
the s channel alone). The main sources of these deviations seem understood. 

As for the internal estimate, we find that for the tightest acollinearity cut both ALIBABA and 
TOPAZO miss terms of the order of 0.2-0.6%, depending on the c.m.s energy. For TOPAZO the 
bulk of the effect is in the non s-channel. From the point of view of TOPAZO the analysis of the 
full Bhabha cross-section is anyhow preferable to the t-channel subtraction procedure, provided that 
the s channel dominates. The size of the uncertainty in ALIBABA resembles the one in TOPAZO, 
although this is somehow accidental in view of the different origins of the contributing effects. The 
ALIBABA program is better suited for the t-channel subtraction procedure, with an uncertainty in 
the non s-channcl components ranging from 1.0 pb to 1.8 pb for the various c.m.s. energies and angular 
acceptances. If subtraction is performed then the TOPAZO program is better suited for the s-channcl 
analysis due to the high precision achieved there. 
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